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Abstract 

The reduced transition probability -B(E2;0^ — > 2f) for 28 S was obtained experimentally using 



Coulomb excitation at 53 MeV/nucleon. The resultant B(E2) value 181(31) e 2 fm 4 is smaller than 
the expectation based on empirical I?(E2) systematics. The double ratio \M n /M p \/(N/Z) of the 
0+ — > 2f transition in 28 S was determined to be 1-9(2) by evaluating the M n value from the known 
B(E2) value of the mirror nucleus 28 Mg, showing the hindrance of proton collectivity relative to 
that of neutrons. These results indicate the emergence of the magic number Z = 16 in the \T Z \ = 2 
nucleus 28 S. 
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Magic numbers characterize the shell structure of fermionic quantum system such as 
atoms, metallic clusters [l| and nuclei ^]. A unique feature of the nuclear system is the 
fact that it comprises two types of fermions, the protons and neutrons, and hence the magic 
numbers appear both for protons and neutrons. Most of the recent studies regarding the 
magic numbers are for neutron-rich nuclei. Disappearance of the conventional magic numbers 
of N=8, 20 and 28 (3 5 1 or the appearance of the new magic number N = 16 has been 
shown. They are associated with nuclear collectivity, which is enhanced, for instance, in the 
neutron-rich N = 20 nucleus 32 Mg caused by disappearance of the magic number [j], flo| . 

The new neutron magic number N = 16 has been confirmed experimentally for 27 Na 



(\T Z \ = 5/2) and more neutron-rich isotones 6l-l8l.llll.ll2|. Its appearance can be theoretically 
interpreted as a result of a large gap between the neutron 0^3/2 and S1/2 orbitals caused by the 
low binding energy 6] and/or the spin-isospin dependent part of the residual nucleon-nucleon 
interaction [l^ . In analogy to the magic number N = 16, the proton magic number Z = 16 
must also exist in proton-rich nuclei. However, it has not been identified experimentally in 
the proton-rich sulfur isotopes. The present Letter reports on a study of the magic number 
Z = 16 at the most proton-rich even-even isotope 28 S with \T Z \ = 2 through a measurement 
of the reduced transition probability Z?(E2;0+ — > 2f). 

The B(E2) value is directly related to the amount of quadrupole collectivity of protons. 
The relative contribution of the proton- and neutron-collectivities can be evaluated using 
the ratio of the neutron transition matrix element to the proton one (the M n /M p ratio) for 
0+ -> 2+ transitions QQ- M p is related to B(E2) by e 2 M 2 = B(E2; 0+ -> 2+). The M n 
value can be deduced from the M p value in the mirror nucleus, where the numbers of protons 
and neutrons are interchanged. If collective motions of protons and neutrons have the same 
amplitudes, the double ratio \M n /M p \/(N/Z) is, therefore, expected to be unity. Deviation 
from \M n /M p \/(N/Z) = 1 corresponds to a proton/neutron dominant excitation and should 
indicate a difference in the motions of p rotons and neutrons. Such a difference appears 
typically for the singly-magic nuclei [lj . [l^ . For proton singly- magic nuclei, the proton 
collectivity is hindered by the magicity, leading to \M n /M p \/(N/Z) > 1. For example, the 
singly- magic nucleus 20 O has a large double ratio of 1.7 ~ 2.2 for the 0+ — > 2\ transition 



17H19]. 



We used Coulomb excitation at an intermediate energy to extract the £>(E2;0+ — > 2^) 
value of the proton-rich nucleus 28 S. Intermediate-energy Coulomb excitation is a powerful 



tool to o 
able 101 



:>tain B(E2) with relatively low intensity beams because a thick target is avail- 



20|. The double ratio \M n /M p \/(N/Z) of the 0+ -> 2\ transition is obtained by 
combining the B(E2) values of 28 S and the mirror nucleus 28 Mg. 

The experiment was performed using the RIBF (Radioactive Isotope Beam Factory) 
accelerator complex operated by RIKEN Nishina Center and Center for Nuclear Study, 
University of Tokyo. A 28 S beam was produced via projectile fragmentation of a 115- 
MeV/nucleon 36 Ar beam from the K = 540 MeV RIKEN Ring Cyclotron incident on a 
531 mg/cm 2 thick Be target. The secondary beam was obtained by the RIKEN Projectile- 



fragment separator (RIPS) 21| using an aluminum energy degrader with a thickness of 
221 mg/cm 2 and a wedge angle of 1.46 mrad placed at the first dispersive focus. The 
momentum acceptance was set to be ± 1%. A RF deflector system 22J was placed at the 
second focal plane of RIPS to purify the 28 S in the beam with intense contaminants (mostly 
of 27 P, 26 Si and 24 Mg) that could not be removed only by the energy loss in the degrader. 
Particle identification for the secondary beam was performed event-by-event by measuring 
time of flight (TOF), energy loss (AE), and the magnetic rigidity of each nucleus. TOF was 
measured by using a RF signal from the cyclotron and a 0.1 mm thick plastic scintillator 
located 103 cm upstream of the third focal plane. AE was obtained by a 0.1 mm thick silicon 
detector placed 117 cm upstream of the third focal plane. The average 28 S beam intensity was 
120 s _1 , which corresponded to approximately 1.9% of the total intensity of the secondary 
beam. The secondary target was a 348 mg/ cm 2 -thick lead sheet which was set at the third 



focal plane. The average 



Three sets of PPACs 



Deam energy at the center of the lead target was 53 MeV/nucleon. 



23j were placed 155.6 cm, 125.6 cm, and 66.2 cm upstream of the 



secondary target, respectively, to obtain the beam trajectory on the secondary target. 



An array of 160 Nal(Tl) scintillator crystals, DALI2 24|, was placed around the target to 
measure de-excitation 7 rays from ejectiles. The measured full energy peak efficiency was 30 
% at 0.662 MeV, in agreement with a Monte-Carlo simulation made by the GEANT4 code, 
and the energy resolution was 9.5 % (FWHM). The full-energy-peak efficiency for 1.5 MeV 
7 rays emitted from the ejectile with the velocity of 0.32c was evaluated to be 16% by the 
Monte-Carlo simulation. 

The scattering angle, energy loss (AE), and total energy (E) of the ejectiles from the 
lead target were obtained by a detector telescope located 62 cm downstream of the target. 
It consisted of four layers of silicon detectors arranged in a 5 x 5 matrix without 4 detectors 
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FIG. 1. Doppler-shift corrected 7-ray energy-spectrum in the Pb( 28 S, 28 S 7)Pb reaction. The fit 
by the response function (dashed curve) and the exponential background (dotted curve) is shown 
by the solid curve. 



at the corners for the first two layers, and a 3 x 3 matrix for the third and fourth layers. 
The silicon detectors in the four layers had an effective area of 50 x 50 mm 2 and a thickness 
of 500, 500, 325, and 500 /im, respectively. The detectors in the first and second layers had 
5-mm-wide strip electrodes on one side to determine the hit position of the ejectiles. The 
AE-E method was employed to identify 28 S. The mass number resolution for sulfur isotopes 
was 0.35 (ler). The angle of the ejectile was obtained from the hit position on the telescope 
and the beam angle and position on the target measured by the PPACs. The scattering 
angle resolution was 0.82 degree. 

The Doppler-shift corrected 7-ray energy-spectrum measured in coincidence with inelas- 
tically scattered 28 S is shown in Fig. [TJ A peak is clearly seen at 1.5 MeV. The spectrum was 
fitted by a detector response obtained by the Monte-Carlo simulation and an exponential 
background. The peak energy was obtained to be 1.497(11) MeV, which was consistent with 



the previous measurement, 1.512(8) MeV, by the two neutron removal reaction on 30 S 25|. 
This peak has been assigned to the transition from the 2^ state to the + ground state 25]. 
In extracting the inelastic cross section, transitions feeding the 2^ state were not accounted 
for, because the proton separation energy of 2.46(3) MeV is relatively low and no higher 
excited states were seen in the present spectrum and the two-neutron removal reaction on 



25| . This was supported by the location of the second excited state in the mirror nucleus 



30g 

28 Mg of 3.86 MeV 



The angular distribution of the scattered 28 S excited to its 1.5 MeV state is shown in 
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FIG. 2. (a) Angular distribution for the Pb( 28 S, 28 S 7)Pb reaction exciting the 1.5 MeV state in 28 S. 
The solid curve represents the best fit with ECIS calculation assuming AL = 2. The dashed and 
dotted curves show the Coulomb and nuclear contributions, respectively, (b) Detection efficiency 
calculated by the Monte-Carlo simulation. 



Fig. EJ^a). Figure Mb) shows the angle- dependence of the detection efficiency for scattered 
28 S obtained by a Monte-Carlo simulation. It took into account the spacial and angular 
distributions of the 28 S beam, the size of the silicon detectors, and effect of multiple scattering 
in the target. The cross section integrated up to 8 degree was obtained to be 99(16) mb 
by taking into account the angle- dependent detection efficiency. The error was nearly all 
attributed to the statistical uncertainties, while the systematic errors of the 7-ray detection 
efficiency and the angle- dependence of the detection efficiency were also included (3%). The 
distribution was fitted by that for an angular momentum transfer of AL = 2, calculated by 



the coupled-channel code ECIS97 26| taking into account the scattering angle resolution. As 
seen in the figure, the AL = 2 distribution well reproduced the experimental one, supporting 
the 2 + assignment for the 1.5 MeV state. The ECIS calculation is almost equivalent to the 
distorted wave born approximation, since higher-order processes are negligible in the present 
experimental conditions. The optical potential parameters were taken from the study of 



the 17 + 208 Pb elastic scattering at 84 MeV/nucleon (22|. The collective deformation 
model was employed to obtain a form factor for nuclear excitation. The Coulomb- and 
nuclear-deformation parameters f3c and (3n were employed to obtain the B(E2) value as 



B(E2) = (3Ze J R 2 /47r) 2 /3£. f3 N is related to (3 C by a Bernstein prescription Q, 

(3 N l + (b F /b F )(M n /M p ) 



1) 



Pc l + {b F /b F ){N/Z) ' 
where b F ,-> is the interaction strength of a probe F with neutrons (protons) in the nucleus. 



b F lb F is estimated to be 0.81 for the inelastic scattering on Pb at around 50 MeV/nuc 



19(| . The M n was deduced from the adopted B(E2) value of the mirror nucleus 28 Mg 



con 



28]. 



The B(E2) value for 8 S was obtained by adjusting /3 C and hence M p with /3n calculated 
by eq. ([T]) to reproduce the experimental angular distribution. The dashed and dotted 
curves in Fig. [2] shows the Coulomb and nuclear contributions, respectively. The use of 



the optical potential determined from the 40 Ar + 208 Pb scattering 29J gave a 5.5% smaller 
B(E2) value. By taking the average of the results with the two optical potentials, the 
£>(E2;0+ —> 2+) value was determined to be 181(31) e 2 fm 4 . The associated error included 
the uncertainty of the measured cross section and the systematic error due to the choice of 
optical potentials. The 5(E2;0+ — > 2±) value for the 24 Mg, a contaminant of the secondary 
beam, was obtained to be 444(66) e 2 fm 4 by the same analysis. This agreed with the adopted 



value of 432(11) e 2 fm 4 



28 



exhibiting the reliability of the present analysis for S. 



The B(E2) and E x (2±) values for Z = 16 isotopes are plotted in Fig. [3]^a) and (b), 



respectively. The filled circles show the present results. The open triangles 
and E x (2^) represent known values for the Z = 16 isotopes up to A = 40 



or B(E2) 



28 



The 



B(E2) value increases from S, the neutron singly-magic nucleus, to S, and decreases 
at 28 S. On the other hand, the 2f energy of 28 S is smaller than those of 30 ~ 36 S. These 
features contradict the empirical systematics. For example, Raman proposed the relation 
B(E2) = (25.7±4.5)£ X (2+V_ 1 Z 2 A- 2 / 3 which is obtained by a global fit to E x (2f) and B(E2) 
in a wide range of nuclei 28[. The shaded band in Fig. [3]^a) represents the B(E2) values 
calculated by this formula. As clearly seen, the present data for 28 S is much smaller than 
the expectation of 472(83) e 2 fm 4 . An explanation of these small B(E2) and E x (2±) is given 
by the hindered proton collectivity and the neutron dominance in the 0+ — > 2f transition. 
A similar mechanism is proposed for 16 C laJQ] and 136 Te Q,Q], where small B(E2) and 
E x (2i) values in comparison with neighboring isotopes are observed. 

Figure Ws) shows the double ratio \M n /M p \/(N/Z) of the Z = 16 isotopes. The filled 
circle and open triangles show the present result and the known values, respectively. They 
are obtained by the B(E2) values of the mirror pairs. The open squares represent the double 
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FIG. 3. Plot of the B(E2;0+ — > 2f) values (a), the excitation energies of 2f states (b), and the 



double ratio \M n /M p \/ 
the USDB interaction 



' N/Z) (c) for sulfur {Z = 16) isotopes. The shell model predictions with 



30 ] are shown by the dotted curves for each quantity. The shaded region 



represents the B(E2) predictions by the empirical B(E2) systematics 
represented by the filled circles. 



281 ] . The present result is 



ratios obtained by the combinations of B(E2) and the result of (p, p') on the nuclei of interest 



361 ] . The ratio for 28 S amounts to 1.9(2) by taking the present result and adopted B(E2) 
of 350(50) e 2 fm 4 for the mirror nucleus 28 Mg 28|. The double ratio of 1.9(2) is significantly 



larger than unity indicating again the hindered proton collectivity relative to neutron and 
the neutron dominance in the 0+ — > 2f transition in 28 S. This hindrance can be understood 
if 28 S is the proton singly-magic nucleus by the Z = 16 magi city. This picture is supported 
by the larger B(E2) value and \M n /M p \/(N/Z) ~ 1 of the neighboring iV = 12 isotones: 



2n 



37|, and 432(11) e 2 fm 4 and 0.95(8) for 24 Mg 



28 



The 



356 e 2 fm 4 and 1.05(6) for 26 Si 
double ratios of 30_36 S are close to unity, as seen in the figure, indicating that the hindrance 
of the proton collectivity does not appear in these nuclei. The large double ratios for 38 > 40 S 



can be explained by the neutron skin effect caused by the Z = 16 sub-shell closure 



36 



The dotted lines in Fig. [3] (a)-(c) show shell model predictions with the USDB effective 



interaction using the effective charges of e p = 1.36 and e n = 0.45 |30|, |40j. The calculation 
shows excellent agreement with the experimental E x (2^) values. The overall tendencies of 
the B(E2) and \M n /M p \/(N/Z) are reasonably reproduced. Especially the sudden decrease 
of B(E2) and increase of \M n /M p \/(N/Z) at 28 S are mostly predicted. It indicates that the 
shell model calculation with the USDB interaction accounts for the phenomena observed 
in the present study. It should be note that the model interprets the N = 16 magicity in 
neutron-rich nuclei with the large Si/2-g?3/2 gap, and hence the Z = 16 magicity in proton-rich 
nuclei is inherent in the model reflecting the isospin symmetry. Slight difference remaining 
between the predictions and the experimental data may require further development of the 
theory. 

In summary, the _B(E2;0+ — > 2+) value for the proton-rich nucleus 28 S was measured 
using Coulomb excitation at 53 MeV/nucleon. The resultant B(E2) value is determined 
to be 181(31) e 2 fm 4 . The double ratio \M n /M p \/{N/Z) for the 0+ ->■ 2+ transition in 28 S 
is obtained to be 1.9(2), by evaluating the M n value from the known B(E2) value of the 
mirror nucleus 28 Mg. These results show a hindered proton collectivity relative to that of 
neutrons in 28 S. It indicate the emergence of Z = 16 magicity in the \T Z \ = 2 nucleus 28 S. 
The systematics of the \M n /M p \/(N/Z) values for the Z = 16 isotopes indicates that the 
hindrance of proton collectivity in proton-rich region appears only at 28 S. 

The authors thank the staff of RIKEN Nishina Center for their work of the beam oper- 
ation during the experiment. One of the authors (Y.T.) is grateful for the support of the 
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